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Abstract 

We review the field theory approach to percolation processes. Specifically, we focus 
on the so-called simple and general epidemic processes that display continuous non- 
equilibrium active to absorbing state phase transitions whose asymptotic features 
are governed respectively by the directed (DP) and dynamic isotropic percolation 
(dip) universality classes. We discuss the construction of a field theory representa- 
tion for these Markovian stochastic processes based on fundamental phenomenologi- 
cal considerations, as well as from a specific microscopic reaction-diffusion model re- 
alization. Subsequently we explain how dynamic renormalization group (RG) meth- 
ods can be applied to obtain the universal properties near the critical point in an 
expansion about the upper critical dimensions dc = 4 (DP) and 6 (dIP). We provide 
a detailed overview of results for critical exponents, scaling functions, crossover phe- 
nomena, finite-size scaling, and also briefly comment on the influence of long-range 
spreading, the presence of a boundary, multispecies generalizations, coupling of the 
order parameter to other conserved modes, and quenched disorder. 
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1 Introduction 



1.1 Percolation Processes 

The investigation of the formation and the stationary properties of random 
structures has been an exciting topic in statistical physics for many years. Since 
it provides an intuitively appealing and transparent model of the irregular 
geometry emerging in disordered systems, percolation has provided a leading 
paradigm for random structures. Bond percolation constitutes perhaps the 
simplest percolation problem. Two fundamental variants of bond percolation 
have been introduced in the past: In the isotropic percolation (IP) problem 
the bonds connecting the sites of a regular lattice (in d spatial dimensions) are 
randomly assigned to be open (with probability p) or blocked (with probability 
1— p), and an agent may traverse an open bond in either direction. In contrast, 
in the directed percolation (DP) problem the open bonds can be passed only 
from one of the two connecting sites, whence the allowed passage direction 
globally defines a preferred direction in space, see FIG. 1. 

In these terms, percolation can be viewed as the passage of an agent through 
an irregularly structured medium, in the sense that the agent can propagate 
through certain regions, whereas it cannot traverse other areas. Though perco- 
lation represents one of the simplest models for random systems, it has in fact 
many applications. Moreover, it yields a prototypical non- equilibrium phase 
transition: For small values of p, regions with open bonds form disconnected 
clusters of typical linear extension ^, wherein the agent becomes localized or 
trapped. For p > pc, a critical threshold value for the frequency of open bonds, 
the probability for the existence of an infinite connected accessible cluster be- 
comes non-zero. Therefore, an agent on an infinite cluster is not confined to 
a finite region of space any more, but may percolate through the entire sys- 
tem. This percolation transition defines a genuine critical phenomenon as en- 
countered in equilibrium statistical mechanics, with ^ providing the divergent 
length scale as Pc is approached. It turns out that many, microscopically quite 
different percolation-type systems share their critical properties either with 




Fig. 1. Directed percolation (DP) in two dimensions: bonds connecting the lattice 
sites can only be passed from left to right. 
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Fig. 2. Two-dimensional critical isotropic (left) and directed (right, with preferred 
direction down) percolation clusters. (Prom Ref. [55]). 

IP or DP. Thus, the phase transitions in IP and DP, the latter distinguished 
from the former through the broken spatial isotropy, define genuine univer- 
sality classes. For reviews on isotropic and directed percolation the reader is 
referred to Refs. [1,2,3]. Critical IP and DP clusters are shown in FIG. 2. 

Broadbent and Hammersley [4] first directed attention to the dynamics of per- 
colation processes. In contrast to random diffusion, in percolation processes 
the spreading of agents is mainly governed by the disordered structure of the 
medium. Consequently, they define stochastic growth processes where agents 
randomly generate offspring at neighboring sites, with rates determined by 
the susceptibility of the medium. The agents themselves decay spontaneously, 
thereby producing debris. In this way the medium becomes exhausted, and 
its susceptibility is randomly diminished. Hence, static DP in c? + 1 spatial 
dimensions can be directly interpreted as a Markov process in d spatial and 
1 time dimension. The preferred (propagation) direction defines the temporal 
axis, and the history of the process yields the debris that forms the {d + 1)- 
dimensional spatially anisotropic directed percolation clusters. In contrast, in 
a dynamical isotropic percolation (dIP) process the agents grow from a seed in 
an expanding stochastic anmihis in d spatial dimensions leaving behind the de- 
bris in their wake. This debris then forms dimensional isotropic percolation 
clusters. 

Remarkably, the renormalized field theory of DP has appeared originally in 
elementary particle physics in the guise of Reg g eon field theory [5,6,7]. Grass- 
berger et al. [8,9] pointed out that Reggeon field theory is not a Hamiltonian 
field theory but constitutes a stochastic process for which they coined the 
name Gribov process. It represents a stochastic version of Schlogl's 'first reac- 
tion' [10]. Subsequently, the formal connection of Reggeon field theory to DP 
was explicitly demonstrated [11,12,13], and Janssen and Grassberger stated 
the DP conjecture: The critical behavior of an order parameter field with 
Markovian stochastic dynamics, decoupled from any other slow variable, that 
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Fig. 3. Directed percolation represented as a forest fire under the influence of a 
storm from left to right; burning trees are depicted as red, and burned dead trunks 
are shown black. 

describes a transition from an active to an inactive, absorbing state (where 
all dynamics ceases) should be in the DP universality class [13,14]. The field 
theory of dIP was initiated by Grassberger's formulation of the general epi- 
demic process (GEP) on a lattice, and his statement that near criticality this 
process produces isotropic percolation clusters [15]. This insight led to the 
construction of the renormalized field theory for the dIP universality class 
[16,17]. 

1.2 Universality Principles of Percolation 

It is tempting to express percolation processes in the language of an epidemic 
disease such as blight in a large orchard, the spreading of bark beetles in a 
forest, or proliferation of a forest fire, as depicted in FIG. 3. Here, the suscep- 
tible individuals or healthy trees form the medium, and the disease, blight, 
beetles, or fire represent the agent. The sick individuals, the befallen or burn- 
ing trees may randomly infect neighboring individuals. The sick individuals 
are allowed to recover, becoming susceptible anew. This behavior defines the 
so-called simple epidemic process (SEP), also known as epidemic with recov- 
ery [18]. Properties of the SEP near the ensuing critical point separating the 
endemic and pandemic phases (whose precise location depends on the suscep- 
tibility of the individuals) shares the universal properties of DP. In contrast, 
in a situation where the sick individuals die out or become immune as opposed 
to susceptible again, the medium is eventually exhausted. This scenario de- 
fines the general epidemic process (GEP), also termed epidemic with removal 
[18,19,20]. Its universal properties are governed by the dIP universality class. 
The statistical properties of the debris clusters that are left behind after the 
disease is extinguished are described by static IP. 

In the recent physical literature, basically two different methods have been 
employed to construct field theories of percolation processes. The first ap- 
proach starts from a specific microscopic stochastic reaction-diffusion model 



4 



formulated in terms of a master equation on a lattice. Basically it consists of 
a description of the random walks generated by diffusing and reacting parti- 
cles. A field theory is then obtained through a representation of the master 
equation by means of a bosonic creation/destruction operator formalism and 
introduction of coherent-state path integrals, followed by a naive continuum 
limit [21,22,23,24,25]. The second method, to be applied in the first part of 
this article, takes a more universal, phenomenological perspective, and can 
be viewed as a dynamical analog of the established Landau-Ginzburg- Wilson 
approach to static critical phenomena. It is based on formulating the fun- 
damental principles shared by all systems in a given universality class, and 
directly employs a mesoscopic stochastic continuum theory. Of course, one 
has to a-priori define the slow dynamic fields in the theory, typically local 
densities of the order parameter and any conserved quantities, as character- 
izing the universality class under consideration. This approach assumes that 
on a mesoscopic scale the motion of the fast microscopic variables can be ab- 
sorbed into the statistical properties of the stochastic equations of motion for 
the relevant slow variables. The corresponding stochastic processes should be 
Markovian provided all slow variables are retained. Taking into account the 
general (symmetry, conservation, etc.) properties that define the universality 
class, the stochastic field theory is constructed using small-density and gradi- 
ent expansions, i.e., long- wavelength and low- frequency approximations that 
focus on the infrared (IR) properties of the system. 

For percolation processes, the fundamental statement regarding the slow vari- 
ables consists in the assertion that percolation near the critical point can be 
described by a Markov process that only depends upon the density n{r,t) of 
activated sick parts of the medium and, in the case of dlP, upon the density 
of the debris m{r,t). In terms of the field n the density of the debris is given 
by m{r,t) = A j'i^n{r,t') dt' , where A denotes an appropriate kinetic coef- 
ficient. Indeed, this assumption that percolation can be described solely in 
terms of n becomes manifest in typical lattice simulations of such processes. 
Near the critical point one observes a patchwork of regions consisting of cor- 
related active neighboring sites and vacuum regions devoid of any activity. 
In stark contrast to this picture, simulations for branching and annihilating 
random walks show anticorrelating behavior, namely widely separated single 
active sites propagating on isolated paths and reacting only after encounter 
[3,26]. 

We now proceed to formulate four principles that allow the explicit construc- 
tion of a mesoscopic stochastic description of percolation processes. We shall 
describe the universal aspects of the dynamics of n beginning with the DP 
universality class: 

(i) The susceptible medium becomes locally infected, depending on the den- 
sity n of neighboring sick individuals. The infected regions recover after 
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a brief time interval. 

(ii) The state with n = is absorbing. This state is equivalent to the extinc- 
tion of the disease. 

(iii) The disease spreads out diffusively via the short-range infection (i) of 
neighboring susceptible regions. 

(iv) Microscopic fast degrees of freedom may be captured as local noise or 
stochastic forces that respect statement (ii), i.e., the noise alone cannot 
regenerate the disease. 

For the dIP universality class, we need to modify principles (i) and (ii) by 

(i') The susceptible medium becomes infected, depending on the densities 
n and m of sick individuals and the debris, respectively. After a brief 
time interval, the sick individuals decay into the immune debris. The 
debris ultimately stops the disease locally by exhausting the supply of 

susceptible regions. 

(ii') The states with n = and any spatial distribution of m are absorbing. 
These states are equivalent to the extinction of the disease. 

It is important to realize that the mechanism (i') introduces memory into 
the stochastic process. Note that we do not explicitly assume diffusive mo- 
tion of the sick individuals. All that is needed is a spreading mechanism for 
the agent, here the disease itself: 9* n(t, r)| spread = / cfirP{r — r')n{t,r') with 
a spreading probability P{r) that is assumed short-range, hence allowing to 
approximate its propagation by simple diffusion. However, this scenario by 
no means excludes react ion- diffusion systems where the "infected" particles 
themselves spread out diffusively, if otherwise the general percolation princi- 
ples are satisfied. 



2 Field Theory Representation of Percolation Processes 

2.1 Dynamic Response Functional for Stochastic Percolation Processes 

We proceed to set up the stochastic equations of motion for percolation pro- 
cesses according to the general principles formulated in the preceding section, 

dtn{v,t) = V{r,t), dtm{r,t) = \nir,t), (1) 
or, in a discretized version (in the Ito sense) {t} {tk} with tk+i — tk + A, 

n{r,tk+i) ^n(r,tk) + AV(r,tk) , m{r,tk+i) ^ m{r,tk) + X An(r,tk) . (2) 

Here, V(r,t) denotes a random density that has yet to be constructed, and the 
second equations for m only appear for the dIP process. For simplicity we shall 
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suppress spatial arguments in the following considerations. The Markovian 
assumption implies that the stochastic variables V{tk) are uncorrelated at 
different times, and that the statistical properties of V(tfc) depend solely upon 
n{tk) and m{tk). As a consequence, the generating function of the cumulants 
of V (i.e., the Laplace transform of the corresponding probability distribution) 
must have the general form 



exp (A^V(ife)n(t,+0) = exp (a^ E ^^.K^^)' ^(^^l) > (3) 

^ k ^ ^ k h=i y- ^ 

where the overbar denotes the statistical average over the fast microscopic 
degrees of freedom. The n(tfc) constitute independent new variables. Of course, 
in the DP process the cumulants Ki are independent of the debris m. 

The statistical properties of the stochastic process are fully encoded in the 
simultaneous probability density of the history {n(to), ^(^i), • • • , ^(^fc), • • •}: 



V{{n{t)}) = ll^[nitk+,) - n{tk) - AV{tk)] (4) 
dn{tk+i) I ^^^^^^^^ [AV(tfe) + n{tk) - n{tk+i] 

V I. 
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Upon inserting (3) and reverting to continuous time, we may then formally 
write 7-" as a path integral over functions h{t): 



CO 



n{n{t)}) = I V[h] exp ldtl^p^Ki[n{t),m{t)]-h{t)dtn{t)y (5) 



This expression is always to be interpreted through the preceding prepoint 
(Ito) discretization; we note that this specifically implies 6{t < 0) = for the 
Heaviside function that enters the dynamic response functions. The integration 
over the 'response variable' h runs along the imaginary axis from — ioo to +ioo. 
However, this path can be analytically deformed in finite regions. Therefore n 
can attain real contributions. 

All cumulants Ki are subject to a fundamental constraint as a consequence of 
the existence of the absorbing state(s): they must vanish for n = 0. Assuming 
that the functional Ki can be expanded in powers of n and m, this implies 
that all Ki{n, m) oc n or higher powers of n. Moreover, as we shall demonstrate 
below, the higher cumulants Ki with / > 3 are irrelevant in the RG sense. Thus, 
the remaining task is to find the functional forms of the mean-field part Ki 
and the Gaussian stochastic force correlator K2. After reducing the statistical 
properties of V to these two terms, the stochastic equation of motion for n 
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may be written in the usual Langevin form (in the Ito interpretation) 



dtn{t) = Ki[n{t),m{t)]+C{t), (6a) 



at)C{t') = K2[n{t),m{t)]6{t-t'). (6b) 

Note that K2 > since our theory is based on a real field n. This is in 
contrast to the path integral representation of the theory of random walks 
subject to pair annihilation which may formally lead to a negative second 
cumulant [24,25]. 

Statistical averages of any functional of the slow variable n can now be per- 
formed by means of path integrals over the fields n and n with the weight 
exp(— JT") [27,28,29]; reinstating the spatial arguments, the integral measure 
here is P[fi, n] = Y[r,t dh{r,t) dn{r,t)/2ni. Retaining only Ki and K2, the 
dynamic response functional J from Eq. (5) becomes 

J = fdt {n{t) \dt n{t) - Ki[n{t),m{t)] \ - I h{tf K2[n{t),m{t)] \ . (7) 
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Here integrations over d-dimensional space are implicit. 

Fundamental principles (ii) and (ii') classify DP and dIP as absorbing state 
systems: Any finite realization will reach an absorbing state in a finite time Tabs 
with probability 1. However, in the active pandemic state above the critical 
percolation threshold p,., T;,bs grows exponentially with the system size 
It is thus appropriate to set Tabs — >■ 00 in the thermodynamic limit L — >■ 
00. The assertion that principles (i)-(iv) lead to the DP universality class is 
known in the literature as the DP conjecture [13,14]. The diffusive spreading 
of the disease renders the cumulants Ki and K2 local functionals of n and 
m, whence a gradient expansion is appropriate. The absorbing state condition 
then inevitably implies 

Xi(n,m) =i?(n,m)n + AV^n + ... , (8) 

where the ellipsis denotes higher-order gradient terms. In the same manner, 

K2{n, m; r - r') = 2 [r(n, m) n + A' (V^n) - A" nV^ + . . .] 5{r - r') . (9) 

The contributions with kinetic coefficients A' and A" can be interpreted as 
purely diffusional noise if A" = 2A'. Yet here these terms simply arise from the 
gradient expansion. We shall demonstrate below that these as well as higher- 
order terms are in fact irrelevant for the critical properties of percolation. 
Naturally, the dependence of m is absent for DP. 

Furthermore, n and m are small quantities near the critical point. This allows 
for a low-density expansion 

R{n,m) ^ -\{t + gin + g2m+ . .}j , r(n, m) = A (g'a . . .) , (10) 
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with g2 = for DP. To ensure stability, the couphng constants gi are assumed 
to be positive; otherwise one has to extend the expansion to the subsequent 
powers. Note that there is no general reason for g^ to vanish as is the case 
for interacting random walks [24,25]. Discarding the irrelevant higher-order 
expansion terms, indicated by the ellipsis in Eqs. (8)-(10), and inserting into 
Eq. (7), we at last arrive at the fundamental dynamic response functional for 
percolation processes: 

J ^ J d'^rdtl^h at + A(r- V^) + >^[gin + g2m - g^nj n-gn|. (11) 

Here we have introduced an additional external source q — q{r, t) for the 
agent. Specifically, a seed inserted at the origin r = at time t = in order to 
initialize a spreading process is modeled by q{r,t) = 6{r) S{t). The evaluation 
of the path integrals is always restricted by the initial and final conditions 
n{t — >• —oo) = = h{t — > +oo), respectively. After eliminating the fast 
microscopic degrees of freedom, and thereby focusing on the IR properties of 
the relevant slow variables, the wave vectors q of the fluctuations of the Fourier 
transformed fields ^q(t) and nq(t) is clearly limited by a momentum cutoff A. 
Correspondingly, those fluctuations are effectively set to zero for |q| A. 



2.2 Master Equation Field Theory for a Specific Reaction- Diffusion Model 



An alternative approach to constructing a field theory representation for a 
stochastic system starts from the classical master equation that defines the 
process microscopically. For reaction-diffusion systems, one may formulate the 
reactions in a straightforward manner in terms of (bosonic) creation and an- 
nihilation operators, and therefrom via coherent-state path integrals proceed 
to a field theory action [21,22,25]. This method invokes no phenomenologi- 
cal assumptions or explicit coarse-graining, but does rely on a the validity of 
a 'naive' continuum limit. One may thus view this procedure in analogy to 
the derivation of the $^ theory from a soft-spin Ising model on a lattice by 
directly taking the limit of zero lattice constant. Universal features do not 
immediately become apparent in such a treatment, whereas the approach laid 
out in the preceding section focuses on general principles and in this sense 
rather corresponds to the Landau- Ginzburg- Wilson model for second-order 
equilibrium phase transitions. The phenomenological construction certainly 
requires considerable a-priori insight on the defining universal features of a 
given system, but its benefit is considerable predictive power. Yet there are 
cases, such as pure pair annihilation processes [23,24] or branching and an- 
nihilating random walks with even offspring [30], for which as yet only the 
master equation formulation has led to successful field theory representations. 
In these processes the random-walk properties of individual particles govern 
the essential physics; hence the microscopic picture is essential and a meso- 
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scopic description simply in terms of the particle density would be erroneous. 
On the other hand, proceeding from the microscopic model via taking a naive 
continuum limit definitely fails, e.g., for the pair contact process with diffusion 
(PCPD) [31,32]. 



We shall consider here a specific reaction-diffusion model that displays an 
active to absorbing state phase transition which, according to our general 
principles, should belong to the percolation universality classes. Yet reaction- 
diffusion systems do of course not provide an exact description of, e.g., directed 
bond percolation or the contact process since they differ manifestly with re- 
spect to the microscopic details. We note that one may also construct a direct 
field theory representation for the pair connectivity of bond-percolating sys- 
tem [12,33]. As usual in reaction-diffusion models, the activated particles both 
represent the diffusing agents and may locally create new offspring. We intro- 
duce the following reaction scheme: 



A A 0, 



K 

A^B, A + B^B, 



(12a) 
(12b) 



supplemented with hopping of the agents A to nearest-neighbors on a d- 
dimensional lattice with sites (i,j, ...) with diffusion constant A. The back- 
reaction in (12a) may also be viewed as effectively capturing mutual exclusion 
on the same site. The reactions in (12b) are specific to dlP where the presence 
of the produced debris B suppresses the agent A. 

The reaction rules (12) are reformulated in terms of a master equation that de- 
scribes the time dependence of the probability P{{n, m}, t) for a given configu- 
ration of site occupation numbers {n} = (. . . , Uj, . . .) and {m} = (..., rrii, . . .) 
of the agent A and the debris B, respectively. The configuration probability is 
encoded in the state vector \P{t)) = I]{„,m}-P({^T', '"T'}, ^) |{^, ''^j) in a bosonic 
Fock space spanned by the basis | {n, m}). These vectors as well as the stochas- 
tic processes in the master equation are then expressed through the action 
of bosonic creation and annihilation operators {a,b} and {a,b}, respectively, 
which arc defined via |. . . , nj, . . .) = |. . . , + 1, . . .), and |. . . , n^, . . .) = 
ni\. . . ,ni — 1, . . .), etc. Subsequently, the master equation can be written in 
the form 



d^\P{t))^-H\P{t)) (13) 
with an appropriate non-Hermitean pseudo- Hamilton operator [21,22,25]. For 



10 



example, the reaction scheme (12) leads to H — H^is + -f^reac with 

<i,j> <i,j> 

j 

+ //(ai -6i)ai + i/(ai - l^Sj^jOj . (14b) 
Here < i,j > denotes a pair of neighboring sites. 

In order to compute statistical averages it is necessary to introduce the pro- 
jection state ('I = (0| njexp(aj + 6j). Using the identity (-j Oj = (-1 = (-1 6^ one 
easily finds the expectation value of an observable A{{n, m}) at time t: 

{A){t)= Y: A{{n,m})P{{n,m},t) = {-\A{{aa,bb})\m), (15) 

{n,m} 

where in the last expression n and m are replaced with the operators aa and 66, 
respectively. The formal solution of the equation of motion (13) reads \P{t)) ~ 
exp(— ii?) I-P(O)). Following standard procedures [23,25] the expectation value 
(15) can be expressed as a path integral 

{A){t) = / Via, 6, a, 6] A{{aa, 66}) exp(— 5'[a, 6, a, 6]) , 



(16) 

with an exponential weight that defines a field theory action. After applying 
a (naive) continuum limit, the action becomes 



S= (frdt 



(a — 1) dttt + A Va • Va + {a — 1) {a — pa + n da) a 



+ (6 - 1) dtb + iji{d-b)a + v{d - 1) 66a 



(17) 



Here, the fields a{r,t), a{r,t), b{r,t), and b(v,t) correspond to the coherent- 
state eigenvalues of the bosonic creation and annihilation operators. The 
integrations are subject to da > 0, bb > 0, and the final conditions are 
a(r, t — > oo) — 1 — 6(r, t — > oo). Thus often a field shift according to a = 1 + d, 
6 = 1 -|- 6 is useful. The new variables 6 then appear only linearly in the ac- 
tion (17). Hence, they can be integrated out, which results in the differential 
equation constraint 

dtb = {p — u db) a , (18) 

and the new action 



S — j d'^r dtd dt — AV^ + {a + ijL — p) + {Ka + ub — pd) + Kda 



a, (19) 



where 6 is given by the solution of Eq. (18). If we now discard all fourth-order 
terms which should become irrelevant after the application of coarse-graining 
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near the critical point, the action (19) attains the same form as the response 
functional (11), with the correct sign of the coupling constants. 



Physically, however, the fields are of different origin: Whereas n in the action 
i7 represents the fluctuating density of the active medium, the field a in (17) 
or (19) results from the bosonic annihilation operators for the random walks; 
therefore, a generally is a complex-valued quantity, and only a has to be real 
and non-negative. Thus, the action (19) derived above really is akin to the 
original Reggeon field theory [5,6,7]. This can be remedied by means of a 
quasi-canonical transformation to proper density variables d — 1+d — exp(fi), 
a — nexp(— n), 6 = 1 + 6 = exp(m), and b — mexp(— m). After integrating 
by parts the action (17) then becomes 



jd'^rdt^i 



n dtu + A Vn • Vn — n{Vn 



+ 



1 — exp(— n) a — pexp{n) + Kn 



n 



+ in dtm + ii 1 — exp(m — n) 



n + v 



1 — exp(— n) mn> . (20) 



Integration over m, followed by the expansion of the exponentials, and dis- 
pensing with fourth-order terms finally leads to 



^' ^ J d'^r dt |n dtU - X n V^n + n{Vh 
-\- n 



. P + CF ^ 

a + iJ, — p) + Kn + u m — n 



n 



(21) 



with dtm = pn. This action again acquires the same form as the response 
functional (11). Aside from the different meaning of the fields, it is reminiscent 
of the action (19) without the fourth-order contributions. However, irrelevant 
diffusional noise arises in the action (21), and the noise term oc comes with 
a slightly different coupling constant that is invariably negative, even in the 
"free" case with p = 0. Hence, as remarked above, the stochastic processes 
described by the actions S and S' are in fact distinct, yet their universal 
features are identical. They will however usually differ with respect to non- 
asymptotic, non-universal details. 



2.3 Mean-Field Theory and Naive Scaling Dimensions 



The first approximation in the evaluation of path integrals generally consists 
of a Gaussian truncation in the action J' with respect to the fiuctuations 
about the maximum of the statistical weight exp(— JT). The extrema are in 
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turn determined from the saddle-point (mean-field) equations 

dt + A(r - V^) + \{gi n + g2m - 2g^ h) 








5J_ 

5h 

Sn 



n 



dt + A(t - V^) + \{2gi n + g2m- gah) 



+ 



\g2 j dt' n{t') h{t') . 



n 



(22a) 



(22b) 



For i — > oo, the stable homogeneous stationary solutions in the case of homo- 
geneous vanishing source g = are n = 0, and n = if r > 0. For r < we 
obtain n = \r\/gi for (72 = (DP), and n = 0. m, = \r\/g2 for gi = (dlP). 
This demonstrates that within the mean-field approximation the critical point 
is located at r = 0, and in its vicinity on the active side the order parameter 
vanishes as n ~ |t|^ with critical exponent (3 — 1. 

Now let us scale spatial distances a; by a convenient mesoscopic length scale 
fj,^^ ^ A^^ Consequently we find the naive scaling dimensions Xt ~ 
T IJ?, and g2/ gi ~ IJ^- Hence, in the asymptotic long-time and large-distance 
limit, T and g'2/s'i constitute relevant parameters, flowing to 00 under succes- 
sive RG scale transformations. The last relation shows that the coupling gi 
becomes irrelevant, provided (72 > (dlP). Since the action J is dimension- 
less, we find from Eq. (11) that hn ~ /x"^. It is characteristic of the field theory 
for spreading phenomena with an absorbing state that the dynamic response 
functional J contains a redundant parameter [34] that must be eliminated by 
a suitable rescaling 



n 



K 



-1 



n 



Ks, m^KS, 



(23) 



with an amplitude K that generally carries nonvanishing scaling dimension. 
Because both the lowest-order non-vanishing coupling constants in the mean- 
field part of the action and the stochastic noise strength are clearly required 
for a meaningful non-trivial perturbation expansion, it is convenient to choose 
K such that the corresponding couplings attain the same scaling dimensions. 
Thus we set 2Kgi = 2K~^g'i = (? in the case of DP, and Kg2 = 2K~^g^ = g for 
dip, where K^^ parametrizes, e.g., the line of DP transitions in the Domany- 
Kinzel automaton [35]. We note that this amplitude characterizes the non- 
universal "lacunarity" property of percolating clusters which tends to zero 
in the case of compact percolation. The dynamic response functional now 
assumes the two distinct forms 



Jdp = Jd'^rdt!^ 
Jdip = Jd'^rdtis 



dt + A(t - V') + y (s - s)]s - A/i 5 }> , 
dt + A(t - V) + ^ {2S -s)]s-Xhs\, 



(24a) 
(24b) 



for DP and dlP, respectively, where \h = K ^q. Hence, a seed at the origin, 
i.e., a sick individual at (r, t) — (0, 0), is represented by A/i(r, t) — K~^S{r)S{t). 
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After fixing the redundancy in this manner the naive scahng of the fields and 
couphngs is uniquely given by 

DP: ~ //(^-'^)/^ (25a) 

dIP: Sr.s^^,M/\ s^pi^'+'y^ g^pt^'-'^'. (25b) 



Prom the scaling dimensions of the couplings g we infer the upper critical 
dimensions dc — 4: for DP, and dc — 6 for dIP. At this point it is also straight- 
forward to show that all higher-order terms in the gradient and density expan- 
sions that we had neglected before in fact acquire negative scaling dimensions 
near dc- This proves that those terms are indeed irrelevant for the asymptotic 
IR scaling behavior. Note that relevance and irrelevance are assigned here 
with reference to the Gaussian theory. After discarding the irrelevant terms, 
the dynamic response functionals (25) display duality invariance with respect 
to time inversion, 

DP: s(t) ^ -s(-t) , dIP: s(t) ^ -S(-t) . (26) 

However, in general these symmetries only hold asymptotically. 



2.4 IR Problems and Renormalization 



The important goal of statistical theories is the determination of correlation 
and response functions (generally called Green's functions) of the dynamical 
variables as functions of their space-time coordinates, as well as of the rele- 
vant control parameters. In a compact form, one attempts to determine the 
cumulant generating functional 



W[H, H\^hij V[s, s] exp - J[s, s] + (//, s) + {H, ~s) 



(27) 



Functional derivatives with respect to the sources H and 



N,N J 



(28) 



H=H=0 



define the Green's functions (here we suppress all space and time coordinates). 
In general, an exact expression for yV[H, H] cannot be found and one must re- 
sort to a perturbational evaluation. Perturbation theory is developed starting 
from the Gaussian contribution exp(— ^0) to the weight, and the subsequent 
expansion of the remainder exp(— ^j), where Ji — J — Jq. 

The different contributions to the perturbation series are graphically orga- 
nized in successive order of closed loops in terms of linked diagrams which, in 
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translationally invariant theories, can be decomposed into one-line irreducible 
amputated Feynman diagrams that represent the building blocks for the vertex 
functions. The generating functional for the vertex functions r[s, s] is related 
to the cumulant generating functional via the Legendre transformation 

r[s, s] + W[H, H] = {H, s) + {H, s) , (29a) 
with s = — , s = ^ , 29b) 

and vice versa. It is easy to see that to zero-loop ('tree') order r[s, s] is identical 
to J'[s, s]. The vertex functions are defined via the functional derivatives 



N] 



= ■■ ^N,N , (30) 

=s=0 



and they are represented by irreducible diagrams if Fi^o = 0. Note that all the 
Fq^at necessarily vanish as a consequence of causality. The Gaussian parts of 
the response functionals Eqs. (24) define the propagator 

(s(r, t)s{r\ t'))o =: G(r -v',t- t') , (31a) 

^^"'^^-4. -io. + A(q2 + r)' ^^^^^ 

where the momentum integral is limited by the ultraviolet (UV) cutoff A. We 
have within the tree approximation 

ri,i(q, oj) — G'i,i(q, —oj)~^ — iu -\- A(q^ -\-t) -\- loop corrections . (32) 

The non-hnear contributions from Ji yield the vertices of the field theory. 

As is common with quantum and statistical field theories, the naive pertur- 
bation expansion, here based on our response functionals (24), is problematic 
for several reasons [36,37,38,39,40,41,42]. To be specific, consider the inverse 
susceptibility x(t, g'. A, s)"-^ = A~^Fi,i(q = 0,a; = 0) = r 0{g'^) that de- 
scribes the linear response to a stationary homogeneous external source. We 
shall be interested in how the structure of the theory depends on the spatial 
dimension d as encoded in the parameter e = ci,. — d. First, the perturbation 
series in in general is divergent. Although it is characterized by a vanish- 
ing radius of convergence, it presumably constitutes an asymptotic series and 
should be rcsummable. Yet more severe problems arise from the massless crit- 
ical limit. The critical point Tc, defined implicitly by x(tc, g'. A, £)~^ = 0, may 
by dimensional arguments be written in the form 

rc^g'/^S{gA-^/^e). (33) 

However, the function S{z, s) is not calculable by means of perturbation the- 
ory, since it diverges in the limit 2; — > for all positive e < 2. The limit 
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S{(},e) =: S{e) does exist for e > 2, and an analytic continuation to e < 2 
defines the dimensionally regularized Symanzik function S{e). But now S{e) 
displays simple IR poles at all £ = 2/n with positive integers n. Note that at 
these discrete points Tc oc gr^". 

Let us now introduce the new variable f = r — Tc. The vertex functions 
r uj} , r) , regarded as functions of f instead of r, arc free of IR singu- 
larities for f 7^ 0, but they require the perturbationally inaccessible function 
(33). In order to eliminate the IR poles completely from the calculation in 
any dimension d < dc, it is possible at this stage to change variables to the 
correlation length squared [39], 



e(f,^,A)^ = ^lnri,i(q,a;) 



(34) 



Since the function ^, which does not display IR poles either, should be a mono- 
tonic function of f, it can be inverted to f = f{C,,g,A) with r{oo,g,A) = 0. 
Finally one may substitute f by f(^, g, A) and thus arrive at vertex functions 
written in terms of ^, i.e., f^^ ^^({q, a;}, ^, gr. A) = ^^({q, a;}, f(^, 5, A), ^, A), 
which are now calculable in perturbation theory and are devoid of IR singu- 
larities for £ > and ^ < 00. 

We return to the consideration of the analytical properties of the response 
function. Dimensional analysis yields 

X(r, g, A)-' = f F(fA-2, ^A^/^ e) , (35) 

with the perturbational expansion 

00 

F{0,z,e)^l + J2fn{e,e)z'^. (36) 

n=l 

For £ > 0, the functions /„ are divergent in the critical limit ^ — > 0, a direct 
consequence of the fact that x^-*^ ~ |f|''' for f ^ with a critical exponent 7 7^ 
1 if (i < dc- Hence, a series in In 9 appears to be produced by the perturbational 
expansion which needs to be properly resummed. 

Wilson's momentum-shell renormalization procedure [43] fully avoids the IR 
problems of the naive perturbation expansion. Furthermore, in contrast to 
renormalized field theory, Wilson's approach does not require the elimination 
of the IR-irrelevant couplings prior to deriving the RG flow equations, and so 
pre-asymptotic critical as well as crossover properties are calculable. However, 
even though Wilson's RG procedure in that sense is superior to the field- 
theoretic method, it is generally not advisable for a systematic calculation 
of universal properties to higher then one-loop order. In contrast, the field- 
theoretic RG method [36,37] does not proceed by successive ehmination of 
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short- wavelength degrees of freedom and a rcscahng of parameters, but instead 
exploits the UV-renormalizability of the perturbation expansion for d < d,. to 
enable a mapping from the critical to a non-critical region in parameter space 
wherein perturbational calculations are unproblematic. Both methods are fully 
equivalent with respect to describing asymptotic and universal features. 

We rename the original bare fields and parameters according to s ^ s, s ^ 
s, T —>■ T, etc. In accord with the symmetries (26) we choose the following 
multiplicative renormalizations 

s = Z^l^s , § = Z^/^s , G£ = Z-^Z^^ZuUn' , (37a) 

X = {ZZ)-^/^Zx\, T = Z^'^ZrT + fc, h^Z^/^Z^^h, (37b) 
Z = Z forDP, Z = Zx for dip, (37c) 

where — r(l -|- £/2)/(47r)'^/^ is a convenient amplitude, u represents the 
dimensionless coupling constant, and r = at the critical point. The renormal- 
ization constants Z... = Z...{u, /j,/ A,e) can be chosen in a UV-renormalizable 
theory in such a way that 

f7v,iv({q,'^},'^,^,A,A) = j^j^^Tff^^{{q,uj},T,u,X,fi) 1 + 0((AC)"'^ 

(38a) 

GAr,iv({r, t}, T, g, A, A) = Z^'/'Z^/'G^^^Hr, t}, r, u, A, /.) [l + 0((Ae)-^ 

(38b) 

with A = 2 -I- 0(£). Within each successive order of the perturbation expan- 
sion, the renormalized vertex functions F^j^ are thereby rendered finite and 
well-defined. Note that the physical bare and the renormalized vertex function 
display the same dependence on the variables (q, u;,r) in the limit (A^) — > cxo 
up to nonuniversal amplitudes. In principle, the infinite cutoff limit is unphys- 
ical and only employed here to develop a systematic RG mapping that works 
effectively to higher orders in the loop expansion. The theory becomes only 
UV renormalizable at the upper critical dimension dc (super-renormalizable 
below dc) because we have previously eliminated all IR-irrelevant couplings 
and shifted the relevant control parameters (here, r) to be zero at the critical 
point. Indeed, the problematic UV and IR singularities are linked precisely 
at the critical dimension d — dc- For d > dc the field theory is IR-finite and 
UV-infinite, whereas conversely for d < dc it is IR-infinite and UV-finite. 

For d < dc, the renormalization factors diverge in three distinct limits The 
first two represent UV divergences, while the third one constitutes an IR 
singularity: In Z...{u, /i/ A, e) — > oo if 

(1) n/A — >• 0, £ = 0, M fixed; 

(2) n/A ^0,e^O,u fixed; 
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(3) ///A — > Iji/h. 0, — > u{l) w*, £ > fixed. 

Here -u* denotes the first nontrivial zero of tlie Gell-Mann-Low function I3{u) 
to be introduced below, and u{l) is tlie solution of the RG flow equation 
ldu{l)/dl = (3{u{l)). Hence, the determination of the Z factors from the UV 
divergences provides us at the same time with important information on the 
critical IR singularities and thereby on critical exponents. This observation 
lies at the heart of the field-theoretic RG method. Explicit calculations of the 
renormalization constants are facilitated if one first takes the continuum limit 
A — > oo with £ > {dimensional regularization) together with the requirement 
that the Z factors absorb just the e poles {minimal subtraction). Notice that 
this minimal dimensional regularization does not at all require an e expansion. 

Yet the continuum limit raises subtle problems in the realm of statistical 
physics [40,42]. For example, one may conclude from the above remarks that 
the full range from zero to infinity of the bare coupling constant g is mapped 
only to the interval [0, m^,] which lies between the UV-stable fixed point m = 
and the IR-stable fixed point u = u^. Values of u larger than are excluded. 
On the other hand, there exist many pre-asymptotic simulational and exper- 
imental results, e.g., for the three-dimensional Ising model near the critical 
point and some polymer systems, which are linked to this strong-coupling re- 
gion that does not include an UV-stable fixed point. In Wilson's RG approach 
this regime can be reached by a choosing suitable initial values of the relevant 
and irrelevant couplings, and a finite UV cutoff A. The UV-renormalized field 
theory should therefore with A — > oo should therefore be accepted as an ef- 
fective field theory where the influences of irrelevant couplings and the flnite 
cutoff are implicitly encoded in the full range of the coupling u also above the 
flxed point. 



2.5 RG Flow Equations and Critical Exponents 

Next, we compare renormalized vertex functions from the same bare theory, 
but which are renormalized at two different external momentum scales and 
Z/i, respectively. The bare vertex functions are held fixed, and we shall hence- 
forth neglect the 0{^{A.^)~^^ corrections. Then it follows from Eq. (38) that 

riv,jv({q' ij) = ^n,n{^^ ""(O) ^ r^,Ar({q> t{i),u{i), x{i), iii) , 

(39a) 

Gn,n{{^^ 0> \ IJ) = ^n,n{u^ ""(O) G'jv,iv({r, t), r{l),u{l), X{1), l/ij , 

(39b) 
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where = X^/^X^/^ and 



x(u{l), uj = lim 



, X(u{l),u) = lim 
V y A— >oc 



Z(M(/),//x/A,i 



A^oo Z{u, jj,/ A,e) 



(40) 

This is the desired mapping, li t/ij, <C 1, the left-hand site of Eq. (39) is 
IR-problematic, whereas for an appropriate I <^ 1 the vertex function on the 
right-hand site is unproblematic if we choose T{l)/{lixf' ~ 1. It is expected 
that \mii^Qu{l) = u^, and the leading critical properties are transformed to 
the functions X and X. 

Casting this mapping in differential form yields the renormalization group 
equation (RGE). We let A — cxo and define the RG functions as logarithmic 
derivatives with respect to the normalization scale /x, holding bare parameters 
fixed, 



k{u) 



du 



dhiT 



51n/x 



an) 



dlnZ 



d\nX 



d In jj, 



d\nZ 



Consider the renormalized Green's functions 



Gn iv({r, t}, T, u, A, n) = lim 

' A— >oo 



(41a) 
(41b) 

(42) 



Z-''/'Z-''/'G^^j^i{r,t},T,g,X,A) 
using the definitions (41) it is straightforward to derive the RGE 

Gj^M^,t},T,u,X,fj,)^0. (43) 



The solution of this partial differential equation is provided by the method of 
characteristics, which solve the ordinary differential equations 



dlnX 

I — — = l{u) , 



I 



dl 
dlnr 

dl 



I 



k{u) , 



d\nX 

dl _ 
din A 



dl 



(44a) 
(44b) 



with the initial conditions X(l) = 1 = -'^(l), — u, t(1) = r, and A(l) = A. 
This yields the general expression 

GAr,iv({r, 0, r, u, A, //) = X{lf/'X{lf/' G^^^{{r, t}, f{l),u{l), X{1), l^i) . 

(45) 

In the critical limit / —>■ the characteristic equations provide us with the 
asymptotic scaling laws. Prom Eqs. (44) we learn that u{l) flows to a stable 
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fixed point u^, given as a zero of the RG beta function f3{u^) — 0, provided 
the first derivative j3'{{u^) > 0. With the definitions 

T] = j{u^) , fj = 7(xi*) , = 2 - k{u^) , z = 2 + C(it*) , (46) 

and employing dimensional scaling, we obtain from Eq. (45) the asymptotic 
form for the Green's functions 

Gjj^l ({r, t}, T, u, X, I,) = A^i^ I'-:- F^,^({///r, FA^X^iH}, r'/'' A^r / fi') . 

(47) 

Here, A, A, Ax, and At- represent four non-universal amplitudes. The sym- 
metries (26) provide us with scaling relations between the critical exponents. 
For DP, Eq. (37c) implies 7 = 7, whereas 7 = 7 + 2C for dIP. Thus we obtain 

fi = r] for DP , 2{z - 2) = r] - ij for dIP . (48) 

The scaling functions Fjv,iv universal, and the exponents 5^ are given by 

Sj^^^ ^ ^{N + N) + zN , p^u ^~'l^^ for dip. (49b) 

Consequently, there are only three independent scaling exponents and v. 
In principle, a similar reduction may be found also for the number of inde- 
pendent non-universal amplitudes. However, the rescaling (23) introduces an 
additional amplitude K related to A and Ahy K"^ — Aj A. 

In summary, we have established that percolation processes near the critical 
point are asymptotically described by universal scaling functions with three 
scaling exponents, but four non-universal amplitudes. The remaining task is 
to explicitly calculate the Z factors and the universal scaling functions in the 
non-critical region by means of perturbation theory. There are two methods to 
avoid the non-calculable Symanzik function. Either one changes variables from 
the parameter r to the correlation length ^, see Eq. (34), which is invariant 
under renormalization: ^(f, g^. A) = ^{t,u,ii), and then calculates perturba- 
tionally the derivative of r by which eliminates Tc [39]. Alternatively, one 
may apply, in addition to the perturbation expansion, a dimensional e expan- 
sion which formally sets Tc to zero, see Eq. (33). 



3 Field Theory of Directed Percolation 

Directed percolation constitutes perhaps the simplest model of a strictly non- 
equilibrium system that displays a continuous phase transition. In this respect 
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t) < (f = G(q,t-t') 

q 

Fig. 4. Elements of the diagrammatic perturbation expansion: propagator (top) and 
three-point vertices (bottom). 

its role is comparable with the paradigmatic Ising model for equilibrium criti- 
cal phenomena. Correspondingly, the DP field theory as given by the response 
functional (24a) can be regarded as the non-equilibrium analog of the (j)^ field 
theory. In the following sections we shall consider the renormalized DP field 
theory and the ensuing asymptotic scaling of important quantities in detail. 



3.1 Perturbation Theory, Renormalization, and Asymptotic Scaling 



As stated before, the perturbation expansion is arranged loop-wise with re- 
spect to the harmonic part of the response functional (24a). In the momentum- 
time representation, the propagator (31b) reads 

G(q, t) = 9{t) exp [-A(f + q'^)t\ , (50) 

where the Heaviside step function is defined with 6{t = 0) = 0. This follows 
from the Ito discretization of the path integral and ensures causality. The 
anharmonic coupling terms in J^p define the elements of the graphical per- 
turbation expansion, depicted in FIG. 4: An arrow marks a s-"leg", and we 
conventionally draw diagrams with the arrowss always directed to the left (i.e., 
we employ ascending time ordering from right to left). The perturbation series 
of the translationally invariant field theory can be analyzed through calcula- 
tion of the vertex functions r^^({q, a;}) that correspond to the one-particle 
irreducible "amputated" graphs. 

Wc are now ready to consider the renormalization of the DP field theory. To 
this end, we evaluate all diagrammatic contributions to the vertex functions to 
a given loop order by means of dimensional regularization, and subsequently 
absorb the UV divergences order by order in the loop expansion into a suitable 
renormalization of the fields and model parameters. We choose the scheme 
(37) and determine the renormalization constants in the minimal subtraction 
prescription. Absorbing the e poles of the naively divergent vertex functions 
into Z factors in fact renormalizes the full theory. The naively divergent vertex 
functions carry non-negative scaling dimensions. In the case of DP, these are 
Ti,! ~ A*^ and ri,2 = — r2,i ~ A*^, schematically illustrated in FIG. 5. 

We will now explicitly determine the renormalizations to one-loop order. The 
primitively divergent one-loop graphs are shown in FIG. 6. We begin by ex- 
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pressing the one-loop contribution to the propagator self-energy, FIG. 6(a), as 
a function of external momentum q and frequency lo: 



6(a) 



2 h iuj/X + 2f + (p - q/2)2 + (p + q/2)2 



2e ^ \2-e 



luj q \ 

^ + +... . 

2A 4 / 



(51) 



Here we have retained only terms linear in lo and q^. These are the contribu- 
tions that display poles in e = 4 — rf. The first term in the brackets has an IR 
singularity at e = 2 (ci = 2). This pole can be eliminated by changing from 
the variable r to the correlation length ^ as an independent parameter. Here, 
for simplicity, we just employ the e expansion. For the extraction of the diver- 
gences of the vertex function ri^2, the external momenta and frequencies may 
be set to zero. The contribution from the diagram in FIG. 6(b) then becomes 



m = 2A^ 



1 



2(f + p2 



(52) 



Combining the zero-loop expressions with the results of this short calculation, 
we obtain the one- loop vertex functions to the desired order in u and q^: 



ri,i = -(l-§^^r-/^)+Aq^(l-§^^r-/^ 



+ Af 1 



and 



e{2-e) 



ri,2 = A^(l-^^2f-/^j + .... (54) 
In order to absorb the e poles into renormalization factors, we employ the 



(53) 






Fig. 5. The naively divergent vertex functions Fi^i and ri^2 
theory. 



-Fa 1 in the DP field 





(a) (b) 

Fig. 6. The one-loop Feynman diagrams for DP. 
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(a) (b) 

Fig. 7. The two-loop diagrams for the propagator self-energy. 

scheme (37), (38a) using Z — Z. To one- loop order we arrive at 



Til = iuj 



'z-^(^ 



e/2 



Z.-^(^ 



u / fl 



s{2-s) \t 



e/2 



8£ \t 

+ ... , 



e/2 



(55) 



and 



2u fnY'"^ 



e \T 



(56) 



Therefore, the e poles are eliminated by the minimal choices Z = 1 + u/Ae, 
Zx = l + u/8e, Zr = l + u/2£, and Z„ = 1 + 2u/s. 

The two-loop graphs for the propagator self-energy and the three-point vertex 
are depicted in FIG. (7) and (8), respectively. The ensuing explicit Z factors 



(a) 



O 

(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



(h) 



(i) 



(k) 



Fig. 8. The two-loop diagrams for the three-point vertex function; the digits indicate 
the number of different possible time orderings. 
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read [13,44,45] 



Z 
Zx 

Zr 

z,. 



^ /13 31 35 

M /16 ^\ ^. o. 

^ 2ii /4 \ 



(57a) 
(57b) 
(57c) 
(57d) 



The RG functions then follow directly from Eqs. (41), 



an)- 



K[U] 



(5{u) 



7 + 10 In- V O(u^) , 

8 V 37 256 ^ ^ ' 

3-u 



4\ 

169 + 106 In- + 0(m=^) 

3/ 128 ^ ' 



u . 



(58a) 
(58b) 
(58c) 
(58d) 



The IR-stable fixed point is determined via /?(«*) = by means of the s 
expansion. 



2e 



1 + 



169 53 , 4\ ^, 2, 

\ In- £ + 0(e^] 

288 144 3/ ^ ' 



(59) 



As a consequence, the general asymptotic scaling resuhs (46) and (47) are 
finally found to carry the ^-expanded special DP critical exponents 

3. 



77 = 


e 

~6 


V28J 


; + - 


z — 


2- 


£ 

12 




'67' 
288 


V — 


1 

2 + 


e 

16 




'107 
.288 


/3 = 


V 

2 


= 1 - 


e ' 
' 6 . 



59 , 4 

144^^3 
17 , 4 



/II 53 , 4\ 

n- £ + £^ 

V288 144 37 ^ 



(60a) 
(60b) 
(60c) 
(60d) 



with = (A^ + ^)/?/i^- 

3.2 Critical Properties of Directed Percolation 

Equipped with the important general asymptotic scaling results for the DP 
Green's functions, which followed directly from the knowledge of the Z factors. 
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we are now in the position to determine the critical properties of dynamic as 
well as static observables. 

3.2.1 Dynamic Observables 

Using the asymptotic scaling form (47), the number N{t,T) of active particles 
generated by a seed at the origin is given by 

N{t,T) = jd''rGi,i{r,t]T,u;\ij) 

= fi-'^A A (q = 0, 1' Ax Xfj,^ t, r^'^Ar rjii^) 

= A^t'^^^^{B^Tt^l''\\). (61) 

Here we have introduced the exponents 

v^^ = zv, es = i/zu , 7' = di/ - 2/? , (62) 

and the non-universal amplitude combinations 

A = A//2A,, Ajv = /x-'^A^AAFi,i(0,l,0), ^ ^x-^A\'''^^ A, . (63) 

^n{x) is a universal scaling function normalized to $Ar(0) = 1. 

As shown in Ref. [46], the survival probability P{t,T) of an active cluster 
emanating from a seed at the origin is asymptotically given by 

P{t, r) = - lim (e-*^^s(r = 0, -t)) = -Go 1 (0, -t;T,k^ 00, u; A, fx) , (64) 

fe— >CXD ' 

with Af — Jd'^xs{x, 0), and where the Green's function Go,i is calculated with 
the response functional J'k = Jj^p + kAf. Note that because of the duality 
invariance (26) the survival probability P{t, r) is in fact related to the mean 
asymptotic density of the active particles p{t, r; po) for a process which starts 
with a homogeneous active density po by 

A P{t,r) ^ A p{t,r; 00) . (65) 

RecaUing the scahng form (47) of Go,i we deduce that 

P(t,T)^Apt-^^^p(BrTt^/''^\), (66) 

with 

6s=f3/ziy, = ^0,1(0, 1,0, cx)). (67) 

Again, $p(a;) represents a universal function normalized to $p(0) = 1. $p(a;) 
tends to zero exponentially for a; — > 00, and asymptotically to Cp for 
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X — oo, where Cp is a universal constant. Thus, we find the percolation 
probability, for r < 0, 

Poo(t) = Aoo|t|^ A^^ApB^^Cp. (68) 



The extension of a active cluster at time t generated by a seed at the origin is 
measured by the radius of gyration R{t, r) of the active particles, as defined 
via 

Jd''rr^G,,,{T,t) ainG'i,i(q,t) 



From the asymptotic scaling law (47) we infer 



(69) 



q=0 



ainG'i,i(q,t) 



_2 d In Fi,i(q, l^A^ Xfi^ t, l-^/^A^ T/^i^) 



q=0 



dq^ 



q=0 



= (Z//)-2 fn{VA^ l-'I^A, rlii^) , (70) 
whence for the radius of gyration we obtain the asymptotic scaling form 

R\t,T) = A^t'^ ^R{B,Tt^l-\\) , (71) 

with 

z, = 2lz, Ap = ii-''A\^fp{\,^), (72) 
where again ^pix) is a universal function normalized to $_r(0) = 1. The four 
quantities A^. Ap, Ap, and Br define a measurable complete set of non- 
universal amplitudes for any system belonging to the DP universality class. 
Once these amplitudes are determined, all other observables attain universal 
values. 

The last dynamic observable we consider here is the active density p{t, r; po) 
for finite initial po, 



p{t, t; Po) = (^1,0(0, t; T, Po, u; A, p) 



(73) 



The initial density po is introduced into the response functional (24a) via a 
source h{r, t) = X~^po ^{t)- No new initial renormalization [47] is involved since 
the perturbation theory is based solely on the causal propagator (the correla- 
tors do not enter), and initial correlations are irrelevant. Hence, according to 
the scheme (37), the renormalization of po is given by 



5n = Z-'/' 



Po 



Po 



(74) 



which leads to an additional derivative term ^'y pod/dpo in the RGE (43). 
Thus a new dependence on X(/)^/^po and Al^^^'^^/'^po/ p'^^'^ arises in Eqs. (45) 
and (47), respectively, which leads to 



p{t, r; Po) = A, t-'^ <^p{BrT t^/^^n , po t'^+'^) , 



(75) 



26 



where $p(a;, y) is universal with oo) = ^p{x) and ^p{x, y) — ^'p{x) y + 

0{y'^). The non- universal amplitudes are 

Ap = {A/ A) Ap, Bp = ii-'^l-'A Al^+'^ , (76) 

and the initial scaling exponent is 

Oi^-nlz. (77) 

3.2.2 Static Observables 

A steady state with DP dynamics can be generated by introducing a ho- 
mogeneous and time-independent external source of activity h in the re- 
sponse functional (24a). In such a steady state one can then measure static, 
time-independent observables, such as the mean density of activated particles 
and their fluctuations [48]. Recall, however, that DP defines genuine non- 
equilibrium systems; hence, there is no fluctuation-dissipation theorem that 
would relate the correlations of the order parameter to its response to an ex- 
ternal conjugated field. Dynamic correlation and response functions therefore 
constitute independent observables. 

In order to obtain the equation of state, i.e., the order parameter M as a 
function of r and h in the steady state, we perform the variable shift s M+s 
and determine M = Gi^q^t, h, u, A, fi) by the 'no-tadpole' requirement (s) = 0. 
After this shift, the linear and harmonic part of the response functional become 



Jdp,o = jd'^rdtS^s 



dt + A(r + gM-\/' 



XgM .2 



+ A 



M{r+'-f)-h]s]. (78) 



As a consequence, aside from the propagator the perturbation expansion is 
now based on a correlator induced by the noise vertex ~ 5^. We will present 
no details of the calculation here, but merely state the equation of state in a 
parametric form, as derived in a two-loop approximation in Ref . [49] : 

T = aR{l-e), M = R>^e, h = bR^ e{2 - 6) + 0{e^) , (79) 

where a and b are non-universal amplitudes that can be related to the previ- 
ously introduced fundamental four amplitudes, and the exponent A is given 

by 

A = /3 + 7, -f = {z-ri)u. (80) 

The parameter R> 0, which measures the distance to the critical point, and 
the parameter range < ^ < 2 are required to describe the entire phase 
diagram in the critical region. 
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As a simple application of this parametric representation we briefly discuss 
the susceptibility x — dM/dh\T, which satisfies a power law for /i ^ 0, 



X = X±\T\ 



-7 



with amplitudes x+ and x_ that correspond to the cases r > and r < 0, 
respectively. We obtain 

^^-^ A^(2-e)+ 2(1-^)2- ^^^^ 

Therefore the universal amplitude ratio X-/x+ can be expressed to order 
in terms of the order parameter exponent (5 as 

^ = 2/? - 1 + 0{s^) = 1 - - [l + 0.067£ + 0(£2)1 . (83) 
X+ 3L J 



The susceptibility x ™ay also be represented by the integral over space and 
time of the Green's function Gi^i: 

X^xjd^'rdt G'i,i(r, t- r, h, u, X, fi) - ItI^^ f±{\r\-^h) . (84) 

In contrast, the mean-square fluctuation %' = {{AN)^) /V of the activated 
particle number are given by 

X' = / dt G2,o(r, t- T, h, u, A, ii) = |T|-^74(|Tr^/i) , (85) 

where the exponent 7' is defined by Eq. (62) and differs from 7. A one-loop 
calculation yields 

with the non-universal parameter 6' and ct = (| + | In 3/2^ £ -|- 0{£^). 



3.2.3 Logarithmic Corrections at the Upper Critical Dimension 

Above the upper critical dimension d^. = 4, where mean-field theory is appli- 
cable, the coupling constant g tends to zero under the renormalization group 
transformation. However, g represents a dangerously irrelevant variable here, 
since it scales various observables, and setting g — rigorously leads either to 
zero or infinity for relevant quantities. The twofold nature of g as both a rele- 
vant scaling variable and an irrelevant loop-expansion generating parameter is 
lucidly exposed by writing the generating functional for the vertex functions 
(29) in the form 

^s,s;T,g]^ g-^<^[gs,gs;r,u]. (87) 
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The expansion of the functional $[<^, 97; u] into a scries with respect to m = 
Gg.^'^g'^ yields the loop expansion. The zcroth-ordcr term g^'^^[gs, gs\{)\ is 
just the response functional (24) itself, i.e., J represents the mean-field con- 
tribution to the dynamic 'free energy' V. The scahng form of the generating 
functional (27) for the cumulants that corresponds to (87) reads 

W[H, H- T, g] = g-''n[gH, gH; r, u] . (88) 

To leading order in the logarithmic corrections, wc may neglect the dependence 
of Q, and $ on u. (For the next-to-leading corrections see Ref . [50] .) 

Solving the characteristic equation (44) for u at the upper critical dimension 
{e — 0) yields to leading order for I — > 0: 

iZ(0 = lM|lnZ|, (89) 

where we introduce the Taylor expansion f{u) — fo + fiU + f2U^ -|- ... for any 
of the RG flow functions / = 7,C,«,/9. Thus, from Eq. (58d), p2 = 3/2 for 
DP. The remaining characteristics are all of the same structure, namely 

dlnO , . . . 

= q(u) . (90) 

Here, Q stands for X, Inf, and In A, respectively, whereas q represents either 
7, or C, as defined in Eq. (41). To leading order, Eq. (90) is solved by 

Q{1) (X llnll-^'/f^^Qil) . (91) 

Combining everything we obtain asymptotically 

W[H,H;T,g;\,fi]:^\lnl\ 



X fl 



lnl\-^/^^H, I In Z I -5/^2^; | In/j-^/V, 0; | lnl\'/'^X, / J , (92) 



where the non-universal amplitudes have been absorbed into the variables. 
Taking the required functional derivatives at ^ = 0, H{r, t) = Xh + po5{t) and 
employing dimensional scaling then yields the asymptotic Green's functions 
with the logarithmic scaling corrections at dc — ^: 

G^,iv({r, 0, u, A, //) oc I In \ In Z| -Vi2)iv+iv 

X F^,jv({^r, f\ lnZ|Vi2i}. r^] lnZ|-VV, hil\-^l^h, T^] InZl-^/^Vo) ■ (93) 
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Prom this general result we find for the dynamical observables 



(Int) 



(lnt)V6/^((lnt)-i/6^^' 

/p((lni)-Vnr), 



N{t,T) 

P{t,r) 

R'{t,T) ^tQnty/'' fp{{lnt)-'/'tT) , 
(lnt)V2 



= {\nty/'pof,{{lnt)-'/HT, {lnt)-'/Hpo) . 
The logarithmic correction for the percolation probability becomes 

Poo{r) = Ce{-T) |r||ln|r||^/3_ 

and the equation of state and the fluctuations in the steady state 
leading order 



X 



A' 



Vr2 + Bh-T 
1- ^ 



In Vr2 + Bh 



1/3 



In Vt2 + Bh 



Vt^ + Bh_ 

where A, A', B, and C are non-universal amplitudes. 



1/6 



(94a) 

(94b) 
(94c) 

(94d) 

(95) 
read to 

(96a) 
(96b) 



3.2.4 Finite-Size Scaling 

We end this section by discussing a method that allows the explicit calculation 
of finite-size effects in DP [51]. We will consider the model defined by the 
action (24a), in a finite cubic geometry of linear size L with periodic boundary 
conditions. Expanding the fields in Pourier modes 



s(r, t)^Yl e'^^'^lq' t) , s(r, i) = ^ e^^-^5(q, t) , 



(97) 



where each component of q only assumes discrete vahics, namely multiples of 
2tt/L, it is clear that the (q = 0)-mode cannot be treated perturbatively at 
the critical point t = h = 0, since the propagator displays an isolated pole at 
q = 0. Therefore, in order to evaluate finite-size effects, one has to construct 
an effective dynamic response functional for the (q = 0)-mode (which subse- 
quently has to be treated non-perturbatively) by tracing out all modes with 
q^O. 

Upon decomposing the fields 

s(r,t) = <l>(t)+v9(r,t), S{r,t) = ^t) + ^{r,t) , (98a) 
with Jd'^rip{r,t) ^ Jd'^rip{r,t) ^0, (98b) 
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one can separate the (q = 0)-modes and from their orthogonal 
complements <f{r,t) and (p{r,t). After performing this decomposition in the 
action J^dp we obtain 



/ 



dt d\ I (p 



(99a) 
(99b) 

Integrating over all modes with q 7^ we arrive at the effective action 

J^? = - In / V[^, ^\ ew{-jS-jSl-jS) = J^i[^, ^ ■ (100) 

The contribution can now be analyzed pcrturbativcly by means of a 

double expansion in powers of the fields $ and $ that arise in the vertices of 

i7j3P , and in the number of loops due to insertion of extra vertices originat- 

(2) ~ 
ing from J/^p. Up to and including terms of third order in and after 

renormalization, one finds [51] 



= L'^ I d J $ [f a* + A f + ^ ($ - I. ) 



(101) 



In dimensions d < A, the functions r(r, L), f{T, L), and g{T, L) display scaling 
properties that follows from the RGE. For c? > 4, the contributions stemming 
from <l>] can be neglected asymptotically, whence f ~ 1, f ~ r, and g — g. 



We proceed by rescaling the fields and the time scale according to 

^{t) ^aM{s) , l>(t) = Q;M(s), Xt ^ ps , 
_ ^-i/2^-d/2 ^ P = 2g-^r^/''L''/\ 

and thereby obtain for the effective response functional 



J^f^ = JdsiM ds + a+ (^M - M 



where the two parameters a and h are given by 

a = 2f- V2 g-^ f L"'^ = ao + a^rh"^ + O ((tL V-)^) , 
h = 2f g-^ hL'^^ (bo + 0{rL^/''))h L^l"" , 

and 



(102a) 
(102b) 



(103) 



(104a) 
(104b) 

(105) 
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For d > 4, we have a oc g~^L'^/'^T + cgl?'^!"^ and h oc g~^L^h. Hence, besides 
the finite-size scahng of the control parameter r and the source h, an additional 
shift oo results from the elimination of the (q 7^ 0)-modes. However, above 
four dimension this shift is proportional to gl?~'^l'^ and represents the leading 
correction. The Green's functions for the spatially averaged variables 
and are now calculated with the response functional [M, M; a, h\ , 
Eq. (103). Consequently 

G^^^m.r.KgM = ([$]^[$]^)(--) = «^+%^({25},a,6) , (106) 

and above four dimension, i.e., in the mean- field region, we obtain asymptot- 
ically 

G^i^m. K g. L) = L-(^+^)'^/^F^^({A^L-'^/2 0, g-'L'l' T, 2g-'L' h) . 

(107) 

The scaling functions F^jv ^® evaluated by an alternative approach. Via 
the equivalence of the response functional (103) with the Ito-Langevin equa- 
tion, which in turn is equivalent to a corresponding Fokker-Planck equation 
for the probability density V{M, s) of the stochastic variable M, we obtain 



In particular, the stationary solution VstiM) for 6 > is easily found. 



M 



(109) 



where C represents a finite normalization constant. All the moments of the 
averaged agent density p can now be expressed in terms of parabolic cylinder 
functions. In particular, for 0? > 4, and at criticality (a = 0), using p oc $ oc 
lj-d/2j^^ one arrives at 

where Ap and Ah are non-universal amplitudes, and T{x) denotes Euler's 
Gamma function. Specifically, Binder's cumulant Q is given by the simple 
expression 

This function has been measured to high precision in recent simulations [52]. 
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4 Field Theory of Dynamic Isotropic Percolation 

We consider now the field theory of dynamic isotropic percolation. We are 
interested in the dynamic as well as static (emerging as t ^ oo) properties 
of isotropic percolation. As mentioned earlier, the static debris clusters after 
the activation process have ceased are described by the usual static percola- 
tion theory near the critical threshold. Of course, we could base our entire 
RG analysis on the full dynamic functional Jdip as given in Eq. (24b), and 
extract the static from the full dynamic behavior by taking the limit i — > oo 
in the end. This would imply, however, that we would have to determine all 
the required renormalizations from the dynamic Fcynman graphs, composed 
of the diagrammatic elements encoded in the action jTdip- Fortunately, there is 
a considerably more economic approach which is based on taking the so-called 
quasi-static limit. Imagine that we initiate a process through external activity 
sources Xh{r,t) = k{r)6{t), localized in time at t = 0; but we are interested 
only at Green's functions of the debris at t = oo. We will see shortly that 
the perturbation expansion simplifies tremendously in this limit. All renor- 
malization factors but one can be calculated directly using this much simpler 
method. Thus, only for the single remaining renormalization do we have to re- 
sort to the full dynamic response functional jTdip. Taking the quasi-static limit 
amounts to switching the fundamental field variable from the agent density to 
the final density of debris ip{r) := S{r, oo) = A /q°^ dts{r,t) that is ultimately 
left behind by the epidemic, and the associated response field (f{r) = s(r,0). 
In particular, the Green's functions corresponding to the correlation functions 
{Hi <^(rj)<^(0)) will turn out to be important for our analysis, since they actu- 
ally encode the static properties of the debris percolation cluster emanating 
from a seed that is localized at the origin at time t = 0. 

4-1 Quasistatic Field Theory 

4.1.1 Quasistatic Hamiltonian 

Following up on the previous remarks, we now proceed to formally take the 
quasi-static limit of the dynamic functional for dIP. The structure of jT^/p 
allows us to directly let 



This procedure leads us from the action (24b) directly to the quasi-static 
Hamiltonian with source k 





(112) 



Hip — 




(113) 
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By just considering the involved time integrations from 'left' (i.e., the largest 
time involved) to 'right' (the smallest time) in any graph, it is easy to see that 
TiiP in fact generates all diagrams contributing to {YliS{ri,oo)YljS{rj,0)) — 
{Hi ¥'('^i) Uj 'Pi^j))- By itself, however, the Hamiltonian (113) is insufficient to 
describe the static properties of isotropic percolation (IP). As a remnant of 
its dynamical origin, Hip must be supplemented with the causality rule that 
forbids closed propagator loops. 

The propagator of the quasi-static theory immediately follows from the Hamil- 
tonian, 



The vertices contained in Eq. (113) are actually identical to those of the DP 
field theory, whence the elements of the perturbation expansion coincide with 
those depicted in FIG. 4 (with the changed propagator (114) and A = 1). 
Correspondingly, determining the renormalization constants of the quasi-static 
theory proceeds as familiar from DP, and is based on the same diagrams, see 
FIGS. 6, 7, and 8. Noticing that our conventions (37) imply that ip = Z^l'^ip 
and (p — Z^^'^'Pi we explicitly find to two-loop order 



These Z factors coincide with those calculated in Ref. [53] for the Potts model 
in the single-state limit. The renormalization constants are known to three- 
loop order [54]. The RG functions appearing in the RGE read to two- loop 
order [53,16] 




(114) 




(115b) 



(115a) 



(115c) 




(116b) 



(116a) 



(116c) 



The IR-stable fixed point, determined as zero of /3(w), is 



2e\ 
= — 




(117) 
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Thus we recover the well-known critical exponents for isotropic percolation: 



2 — k(u^ 



e 



1 



206 



21 

1 _ 1 

~ 2 ^ 84 

e 



441 

he 



e + 0{e'') 



. 589 ^, 2^ 



2205 



6 



61 



(118a) 
(118b) 
(118c) 



Ji..l.2 Static observables 

Let V{S)dS be the measure for the probability that the cluster mass of the 
debris generated by a seed at the origin is between S and S + dS. We obtain 

V{S) = (^6(^Jd''r^{r)-SyMm]) ■ (US) 

For the probability density V{S) of large clusters with S ^ 1 we may use the 
expansion of the exponential to first order (higher orders asymptotically only 
lead to subleading corrections) and find 

V^{S) = (^s(^ld\-ipir)-S^ipiO)) ■ (120) 



The percolation probability P^o is defined as the probability for the existence 
of an infinite cluster generated from a single seed. Hence 



Poo = l- lim / dSe-'''^P{S) 

c^0+ JO 



1 — lim ( exp 

c^0+ \ 



(^(0) -c Jd'^rip{r) 

Via expanding exp[<^(0)] we arrive at the asymptotic form [46] 

~ - lira (^(0) e-^^) , 



^0+ 



(121) 



(122) 



where Ai = J d'^r (p{r). The virtue of this formula is that it relates the percola- 
tion probability in an unambiguous manner to an expression accessible by field 
theory. For actual calculations the term exp(— cM) needs to be incorporated 
into the quasi-static Hamiltonian; i.e., we must replace the original Hip with 



(123) 



Here, c(x) = c plays the role of a source conjugate to the field ip. Whereas in 
general {(p) = by causality if c = 0, the limit c — > 0+ leads to a non- vanishing 
order parameter P^o in the active phase with spontaneously broken symmetry. 
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In terms of averages (• • • )c with respect to the new Hamiltonian (123), we may 
now write 

Poo = - lini (<^(0))e = -^0,1(0; r, c ^ 0+, u, /x) . (124) 

c^0+ 

With the aid of Eq. (124) and the scahng form (47) we readily obtain that 

Poo-^e{-T)\T\^. (125) 

In order to examine the scahng behavior of P(S'), we consider its moments 

IS'') = / dSS^V(S). (126) 
Using Eq. (120), our general scaling result (47) implies 

{S'^) ~ I {d'rf G,,i({r}, 0, r) ~ \t\P-^^'^^-P) , (127) 
which tells us that 

V^{S, t) = Snsir) = S^'^- /{rS'^") , (128) 

where ^5 is the number of clusters of size S per lattice site. These cluster 
numbers ns play an important role in percolation theory, which conventionally 
employs the scaling exponents Tp and ap as defined in Eq. (128). In terms of 
our earlier critical exponents, these are 

It follows from Eq. (127) that the mean cluster mass (S) (taken as average of 
the finite clusters) scales as 

= M(t) =Mo|t|-^, 7 = di/-2/3. (130) 

Next we consider correlation functions restricted to clusters of given mass S. 
In terms of the conventional unrestricted averages with respect to Hip, these 
restricted correlation functions can be expressed for large ^S" as 

Ci^\{r}, r) = ((^(rO . . . ^{v^) 8 (^J d'r (^(r) - s) (^(0))^'°"^ . (131) 
Their scaling from can again be read off from Eq. (47), 

cSf)({r},T) = |T|^^+'^^Fjv({|Trr}, \t\''''-^S) . (132) 
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By means of these restricted Green functions we can write the radius of gyra- 
tion (i.e., the mean-square cluster radius) of clusters of size S as 



2dJddrCf\r,T) 



whence with Eq. (132) 



Rl^S^/'^ffRirS'^-), (134) 

with the fractal dimension 

Df^d-(3/u. (135) 



We conclude this section by considering the scaling behavior of the debris 
statistics if the initial state is prepared with a homogeneous seed density po = 
h. As discussed above, at the level of the quasistatic Hamiltonian Tiip such 
an initial state translates to a further additive contribution —h Jd'^rip{r). Our 
general scaling form (47) implies that the correlation functions of the densities 
</?(r) for the case of a homogeneous initial condition behave as 

Gj^Mi^}, r,h)^j:^ jid'ff G^^^^Jfr}, {?}, r) 

^ |^|/3(iv+iv)^±|^||^|.^|^ |r|^-<^^/i) . (136) 

It is obvious that the initial seed density h plays the role of an ordering field. 

It corresponds to the ghost field of conventional percolation theory [1] . Hence, 
the Green functions do not display critical singularities as long as /i > 0. For 
a homogeneous initial condition the appropriate order parameter is given by 
the debris density 

M = {ip{v))n = Gi,o(0,T,/i) = \T\^ff{\Tf-'"'h) . (137) 

In the non-percolating phase (r > 0) , the order parameter M is linear in h for 
small seed density, with a susceptibility coefficient that diverges as r — > 0, 

M(t > 0, /i) - r-^/i , (138) 

with the susceptibihty exponent 7 given in Eq. (130). At criticahty (r = 0) 
the order parameter M tends to zero for /i — > as 

M(r = 0, /i) ~ (139) 

with the exponent 

.= ^ = 1.1. (HO, 
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Finally, in the percolating phase (r < 0) the order parameter becomes inde- 
pendent of the initial seed density in the limit h ^ and tends to zero with 
T according to 



M(t < 0) 



~ T 



(141) 



which of course is just Eq. (130). Eqs. (125) and (141) establish explicitly 
that the two distinct order parameters, namely the debris density M and the 
percolation probability Pqo, both scale with the identical exponent /3 near 
criticality. 



4- 1-3 Crossovers 

At this stage we are in a position to discuss crossovers between the two fun- 
damental percolation transitions [55,56]. We have estabhshed earlier that the 
introduction of the memory term oc nmn in the response functional (11) 
changes scaling behavior from DP to dIP. Hence, memory is clearly a rele- 
vant perturbation (in the RG sense) in DP and the process described by Jdp 
asymptotically crosses over to the process described by jTaip, both in d spatial 
and the single temporal dimension. We have also argued above that the static 
asymptotic behavior of isotropic percolation in d dimensions is described by 
the quasi-static Hamiltonian Tiip. If we now imagine a forest fire as a spe- 
cial realization of isotropic percolation and introduce a strong wind blowing 
in a preferred direction e, such a directional disturbance may be described 
by an additional term oc (f{e ■ V)(/7 in Eq. (113). Upon comparing with the 
longitudinal part (^(e ■ V)^</? of the Laplacian contribution, we see that this 
additional directional term is more relevant for the IR behavior. That is, if we 
scale the transverse part of the Laplacian as V^^ = — (e • V)^ ~ /i^, the IR- 
scahng dimensions become (e • V) ~ /x^ and (e • V)^ ~ /x^. As a consequence, 
the contribution (p{e ■ V)^</? is asymptotically irrelevant, and the quasistatic 
Hamiltonian becomes 



where c is a new parameter that describes the directionality of the percolation 
in d dimension. It is now straightforward to show by means of simple rescal- 
ing that by identifying the longitudinal spatial direction with 'time' and the 

{d — 1) transverse subspace directions with a new space, Tijp transforms into 
the response functional for DP, namely JTdp, in {d — 1) spatial and 1 tem- 
poral dimension. Thus, we have established the following schematic crossover 
scenarios: 




(142) 
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DP dip IP DP 

quasistatic 

^DP memory J^ip "^ip directionality Jr>p 
^ limit ^ 



d+1 d+1 d (d-l) + l 

We note that one has to apply some care when trying to capture these crossover 
scenarios by means of the dynamic RG, since both the transitions from DP 
to dip as well as from IP to DP involve different upper critical dimensions. 
Therefore a mere e expansion about the dc of either theory cannot possibly 
access the opposite scaling limit. One can, however, work out the RG flow 
functions that incorporate the entire crossover regimes at fixed dimension d, 
provided one employs the correlation length ^ as independent variable rather 
than T in order to eliminate IR singularities. For more details on this method to 
describe the crossover from isotropic to directed percolation, see Refs. [55,56]. 



4.2 Dynamic Observables 



In order to investigate scaling properties of genuinely dynamical observables 
of dIP, one must resort to the full response functional (24b). For the deter- 
mination of the last independent renormalization factor Z to two-loop order, 
one needs to evaluate self-energy diagrams such as FIG. 6(a) and FIG. 7, 
but now with one temporally delocalized leg of the vertex corresponding to 
the coupling oc sSs in J'dip- We thus obtain from the renormalization of the 
derivative dTi^i/duj\ 



q=a)=0 


[16] 
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(ZZ f^ = 1 + I + - ^ + 5 ln4 - 9 In 3) ^ + Oiu^) . (143) 

In combination with Eqs. (115a) and (37), Eq. (41) then yields the additional 
RG function 

7ii / 1 747 \ 7/^ 

^^""^ ^ "12 + + 91n3 - 51n4j- + 0{u') . (144) 

Inserting the fixed point value (117) into Eq. (46), we find for the dynamic 
exponent 



„ e [937 9 ^ -2 



6 



+ — (51n4-91n3) 



588 98 



'^ + 0(s'). (145) 



The scaling form of the survival probability is equal to the expression (66), with 
the spreading exponent 5s = and the longitudinal correlation exponent 
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z/|| = vz that follow from the fundamental exponents listed in Eqs. (118) and 
(145). Likewise, the radius of gyration is given by Eq. (71) with the spreading 
exponent 



2 . £ 



1231 9 



294 49 



+ — (51n4-91n3) 



-.2 



^ + 0(5^). (146) 



Remarkably, this two-loop approximation for Zg differs only by 1-2% from 
simulation results obtained in the physical dimensions d — 1 and d — 3 [1,17]. 

The number of active particles generated by a seed becomes 

N{t,T) = t'^-Viv(Ti'/''ii) , ds = {dv - 2(3)/v\\ . (147) 

The active density pit, r, po) at time t, initialized by a finite homogeneous 
density po, is 

p(t, T, h) = t-^'-^fp{T t^/^ii , h t'^i+^^+i) (148) 
with an initial scaling exponent 

^, = (2-^-77^)/^. (149) 



5 Conclusions and Outlook: Other Classes of Percolation Processes 



In this overview, we have studied the field theory approach to percolating 
systems. Based on the fundamental and universal features of the simple and 
general epidemic processes, we have constructed a mesoscopic description in 
terms of stochastic equations of motion, which we subsequently represented 
through a path integral with the dynamic response functional serving as the 
appropriate effective action. We have also commented on a more microscopic 
representation that starts from the classical master equation of a specific real- 
ization of such processes. In the bulk of this paper, we have provided a detailed 
description of the analysis of the ensuing stochastic field theories, from ba- 
sic scaling properties to the perturbation expansion and UV renormalization, 
and explained why and how one may therefrom infer the correct asymptotic 
IR scaling behavior. We have derived a number of scaling relations for the 
critical exponents of directed (DP) and dynamic isotropic percolation (dIP), 
and explicitly demonstrated how these are linked to the large-scale, long-time 
properties of numerous static and dynamic observables. For the case of dIP, 
we have also derived the effective quasistatic field theory, which yields the 
scaling exponents of isotropic percolation. We remark that non-perturbative 
RG methods have recently been applied to the DP field theory as well [57]. 

Naturally, there are various possible extensions of the above models, some of 
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which lead to novel critical properties. We end this review with a brief discus- 
sion of some interesting modifications of the standard percolation processes. 



5.1 Long-Range Percolation 



In the standard version of the percolation processes in the language of an 
infectious disease, the susceptible individuals can become contaminated by 
already sick neighboring individuals. At the same time sick individuals are 
subject to spontaneous healing or immunization. In more realistic situations, 
however, the infection could be also long-ranged. As an example, envision the 
spreading of a disease in an orchard where flying parasites contaminate the 
trees practically instantaneous in a widespread manner, provided the time 
scale of the parasites' flights is much shorter than the mesoscopic time scale 
of the epidemic process itself. Thus following a suggestion by MoUison [19], 
Grassberger [58] introduced a variation of the epidemic processes with infec- 
tion probability distributions -P(-R) oc l/R^^" which decay with the distance 
R according to a power law. We will refer to such long-range distributions as 
Levy flights., although a true Levy flight is deflned by the Fourier transform, 
P(q) oc exp(— felql'^) [59], and only Levy exponents in the interval < cr < 2 
give rise to positive distributions [60] . 

The spreading probability in this situation is rendered non-local, 

dn{r, t) 



dt 



spread 



= j(fr' P{\r-T'\)n{Y\t). (150) 



After Fourier transformation of this equation, and after applying a small mo- 
mentum expansion, we arrive at 



an(q, t) 



dt 



spread 



Po-P2q^+Paq'' + 0(q\q") n(q,t), (151) 



where the analytical terms stem from the short-range part of P{R), whereas 
the non-analytical contributions arise from the power- law decay. The constant 
Po is included in the reaction rate as a negative ("birth") contribution to r, 
while p2 q^ represents a diffusional term. In order to decide which of the terms 
in Eq. (151) are relevant, one has to compare with the scaling behavior of 
the Fourier-transformed susceptibility x(q, a;) oc where f] denotes the 

anomalous fleld dimension within the short-range spreading theories deflned 
by the response functionals (24), i.e., f] is given by Eqs. (60a) and Eq. (118a), 
respectively. If o" < 2 — the parameter constitutes a relevant perturbation 
and must be included in a renormalization group procedure. Upon taking the 
leading non- analytical term in Eq. (151) into account, the harmonic part of 
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the response functionals (24) changes to 



Jo[s,s]^ Jd'^rdtl^s(^dt + X T-V^ + f(-V^y^^ yy (152a) 



The mathematical treatment of this field theory with both gradient terms is 
somewhat delicate. Therefore, DP and dIP with long-range spreading were 
originally studied via Wilson's momentum shell RG method, since this ap- 
proach can handle relevant and irrelevant contributions on equal footing [61]. 
Within a one-loop calculation, it was shown that the spreading is dominated 
by the Levy flights for a < 2 — fj, and the scaling exponents change con- 
tinuously to their short-range counterparts aX, a = 2 — fj. The application of 
renormahzed field theory has to distinguish between two separate cases. If 
2 — a — 0{e), one must apply a double expansion in e and a — 2 — a, follow- 
ing the work by Honkonen and Nalimov [62]. Thereby the renormalizations 
were obtained to two-loop order [63], and the one- loop results indeed corrob- 
orate the findings within the momentum shell approach with respect to the 
crossover to short-range percolation. 

In the case 2 — a — 0(1), the diffusional term becomes IR-irrelevant, and 
must be removed in order to obtain a UV-renormalizablc field theory in the 
infinite-cutoff limit. Then, by rescaling of time, / can be set to 1. The usual 
scaling r ~ yields g'^ ~ /x^, and (At)^^ ~ r ~ /x*^, where s ~ s ~ jjfi/'^ with 
£ = 2(7 - d for DP and s ~ n^'^-^^'^, s ~ with £ = 3(j - for dIP. Thus 

we infer dc = 2a and d^ = 3a, respectively, to be the upper critical dimensions. 
The propagator is now G(q, t) = 6(t) exp[— A(r + q'^)t]. The vertex functions 
are analytical functions of the external momenta and frequencies as long as 
T > 0. Thus, the non-analytic Levy flight term in Eq. (152a) proportional to 
(_V2)^/2 does not require renormalization, whence besides Z = Z as usual 
we find for DP that Zx = 1 exactly, while for long-range dIP Z = Zx — 1. 
It turns out that counterterms are only needed for vertex functions with zero 
external momenta. The following identity 

/ du fian = I ^''-'^'^ ^ / d^'k m , (153) 

where d' = 2d/a, is useful for the explicit computation of the Z factors, from 
which subsequently the RG functions (41) are found. To one-loop order for 
long-range DP, the result is 



-€+ — + 0{U^) 

4 



u , (154a) 



k{u)^'^ + 0{u'), ^u)^^{u)^au)^-l + 0{u'), (154b) 
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whereas for long-range dIP 



u . 



(155a) 



3u 



k{u) ^u + 0{u^) , 7H=0, ^{u)^2C{u)^-— + 0{u^). (155b) 



At the long-range DP and dIP fixed points — Ae/T + 0{e'^) and 
e/4: + 0(£^), respectively, we then get the critical exponents 



long-range DP: fj^rj^z-a^-- + 0{e ) , - = cr - y 0{e ) , 



(156a) 
(156b) 



3£ Is 
long-range dip : rj ^ 2(z - a) ^ -— + 0(s^) , - = a-- + 0(£^), 

8 4 



/3 = .^ = l-^ + 0(£^), r) = 



(156c) 
(156d) 



5.^ Percolation Boundary Critical Behavior 



Within the field theory formulation, one can also investigate the infiuence 
of a spatial boundary on critical behavior [64,65]. Generally, for percolation 
processes four possible scenarios can be envisioned: The boundary remains 
inactive, whereas the bulk is critical {ordinary transition), the boundary is 
active, the bulk is critical {extraordinary transition), the boundary is critical, 
but the bulk inactive {surface transition), or both boundary and bulk are crit- 
ical, which represents a multicritical point (the special transition) [66,67,68]. 
Let us consider a semi-infinite geometry, with a boundary plane at z = 0; in 
this situation wc need to supplement the dynamic response functionals (24) 
with the surface action 

^surf = jd'^~^rdtXTss{z = Q)s{z = 0), (157) 

and impose the boundary condition dzs\z=o — Tg s{z = 0) . Naive power count- 
ing yields that the new parameter is relevant, whence its RG fixed points 
are either or ±oo. The ordinary transition scenario corresponds to +oo 
and the special transition to — 0. The presence of the boundary implies 
a different scaling behavior of the fields near the surface as compared to the 
bulk. For example, the surface order parameter acquires new critical exponents 
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DP: = + = + (158a) 

dIP: M°^ = ^-^ + 0(.^), /3f^ = l-| + 0(e^). (158b) 

5.3 Multispecies Directed Percolation Processes 

One may readily generalize the previous mesoscopic description of the simple 
epidemic process to multiple activity carriers a = 1, 2, ... in order to capture, 
say, a variety of population species near an extinction threshold. In the spirit 
of Sec. 2.1 one thus arrives at the coupled DP Langevin equations [44] 

dt Hoc = K V^^a + -Ra[{ria}] + Ca , (159a) 
Ra[{na}] = -Aa (ja + - X! 5'a/3?^/3 + •••), (159b) 

with the stochastic noise correlations 



Ca{t) Ut') = Xa Qa n^{t) 6{t - t') . (160) 

Although this coupled multispecies systems appears to be very rich, it turns 
out that in fact all the ensuing renormalizations are given precisely by those 
of the single-species process, and hence the critical behavior at the extinction 
threshold is quite remarkably just that of DP again [44]. In addition, this 
generically universal model displays an instability that asymptotically leads 
to unidirectionality in the interspecies couplings. A special situation arises 
when several control parameters Tq, vanish simultaneously, which implies mul- 
ticritical behavior of unidirectionally coupled DP processes [69]. In the active 
phase, one then finds a hierarchy of order parameter exponents Pa with 

/3i=/3dp = 1-^ + 0(£2), p, = ^-^ + 0{e'),..., /5, = 1 - 0(e) . 

(161) 

In addition, one can show that the crossover exponent associated with the 
multicritical point is $ = 1 to all orders in perturbation theory [44] . There re- 
mains, however, an unresolved technical issue originating from the emergence 
of a relevant coupling that enters the perturbation expansion [69] . Quite anal- 
ogous features also characterize multispecies generalizations of the general 
epidemic process or coupled dIP processes. 
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5-4 Directed Percolation with a Diffusive Conserved Field 



Several years ago Kree, Schaub, and Schmittmann [70] introduced a model 
that consists of the two-species reaction-difTusion system S — > 2S, 25 — > 
and S + C — > C with unequal diffusion constants species B and C. The active 
particles B, whose density we set proportional to s(r, t), become poisoned by a 
diffusing conserved quantity C with density c(r, t). In a mesoscopic description, 
this KSS-model is described by the coupled Langevin equations 

X-^dtS^W^s- (t + ^s + fey + Cs, (162a) 
7-'9tC = V^c + Cc, (162b) 

with positive parameters A, 7, and g. The stochastic forces (i{r, t) must respect 
the absorbing state condition as well as the conservation property. Hence their 
Gaussian fluctuations are given by 

C(r,t)Cs(r^tO = \-'~g s{r, t) S{r - r') 5(t - t') , (163a) 
Ur,t)UT',t') = r\-^') S{r - r') S{t - t') , (163b) 
Cs{r,t)Cc{r',t')^0. (163c) 

More recently, van Wijland, Oerding, and Hilhorst studied the two-species 
reaction-diffusion system A + B ^ 2B and B ^ A with unequal diffusivities 
[71]. Here, the total number of A and B particles constitutes a conserved 
quantity C. Elimination of the A density in favor of the density of C then 
recovers the stochastic equations of motion (162), (163), yet now with a cross- 
diffusion term in Eq. (162b): 

^-'dtC^V^{c-as) + Cc: (164) 

where (7 is proportional to the difference of the diffusion constants of the B 
and A species. 

Both the KSS and the WOH model represent generalizations of DP via includ- 
ing the influence of a non-critical conserved quantity in the critical dynamics 
of the agent, akin to the generalization of the relaxational model A of near- 
equilibrium critical dynamics to model C [72]. Therefore we propose the label 
DP-C for this modification of the DP universality class. Following the same 
lines that resulted in the response functional (11), and omitting IR-irrelevant 
terms, we obtain the renormalizable response functional of the DP-C class 
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corresponding to the above Langevin equations, 
J = Jd'^rdtls 



dt + A(r - + /c) + ^ {gs - ~gs) 



+ c 



dtc- -fV'^{c- as) 



- 7 (Vc)^ . (165) 



Our first observation concerns the stabihty of DP-C. Consider the mean-field 
approximation for the stationary state. Eq. (164) imphes that c(r) = as{r)+Co- 
Including the constant cq in the variable r, we arrive at 

V's- (r+^s)s = 0, g,s^9 + 2af, (166) 

which demonstrates that the stable homogeneous solution s = for posi- 
tive r and g^e becomes unstable for g^s < 0. Higher orders of s should then 
be included in the density expansion, and one would expect the transition 
to become dicontinuous or first-order. A continuous second-order transition 
therefore requires that the constraint g > —2af be satisfied. An other qual- 
itative view on this instability is illustrated by the following consideration. 
Assume / > in the following, i.e., the density c operates as an inhibitor. If 
now an enhancement of s is created by a fiuctuation in some region of space 
the current contribution jcross — 7cVs shows that the inhibitor fiows into this 
region if a > and subsequently reduces the fiuctuation. However, if o" < 0, 
the inhibitor is reduced by the fiow out of this region, and the fiuctuation of 
s becomes increasingly enhanced. 

Rescaling the fields c, c, and the parameter cr, we may set 7 = 7. The response 
functional (165) possesses the following symmetries under three transforma- 
tions that involve constant continuous parameters ai. 



I 

II 
III 



c — > c -|- CKi ; (167a) 
c— >c-|-q;2, T^T — fa2] (167b) 
s^ags, s^ag^s, (7 ^ V , g^a^^g, g^a^g. (167c) 



Moreover, J' is invariant under the inversion 

IV: c^-c, c^-c, a^-a, f ^ -f . (168) 
In the particular case a — 0, the time refiection 



V: ^g/~gsir,t)^-^~g/gsiv,-t), (169) 
c(r,t) ^c(r,-t), g(r,t) ^c(r,-t) -g(r,-t) (170) 

yields an additional discrete symmetry transformation. The invariance with 
respect to the symmetry V distinguishes the special KSS model from the 
general DP-C class. 
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Symmetry I results from the conservation property of the field c. Symmetries 
III and IV show that dimensionless invariant coupling constants and param- 
eters are defined by u = gg fi^"^, v = w = agf , and the ratio of 
the kinetic coefficients p = 7/A with e — A — d. Symmetry III imphes that 
the response functional J again contains a redundant parameter which can 
be fixed in different ways. Dimensional analysis and the scaling symmetry III 
applied to the Green functions Gjv,iV;M,M — ([s]-'^[s]-^[c]-^[c]^) gives 

Gn,n;m,m = "3'"^G'jv,jv;M,M({r, t}, T, Qig V, ^3 §, a^^g, f, A, 7, //) (171a) 

= (7^"^FjV,iV;M,M({/^I", ^HH}, U, V, W, p) (171b) 

^{g/gf^ ^^^^^^,7V;M,m({/^'^'7/^^0>/^"Vm,^,«^,P) , (171c) 

wherein the UV singularities and critical properties reside in the functions 
Pn,N;M,m and ^N,N;M,M' respectively. If o" = 0, only Eq. (171a) can be used, 
and it is natural to apply a rescaling that leads to g = g, which is fixed also 
under renormalization by the time inversion symmetry V. However, this is not 
the case in the general situation a ^ 0, whence it is more natural to absorb 
(7 into the fields s and s. In other words, one may then hold a — a constant 
under renormalization. 

It is easily seen that loop diagrams do not contribute to the vertex functions 
N-M M with M > 1, whence we infer c — c, c — c, 'j — ^, and as — as. We 
define the remaining renormalizations via 

s = Z^/^s, l = Z^/^S, X = Z-^^^Z-^/^ZxX, r = Z^^ZrT + Tc, (172a) 
G^lg = Z-'/^Z-'/'Zf (Z^u + Ywy , gJ^ = Z^^Z^vii' , (172b) 
G.alf = Z-V2z-V2^-2^^^^^ ^ zl = Z^Z, , (172c) 

where the last setting can be implemented through an appropriate choice for Y . 
For 0" ^ 0, the non-renormalization oi s {Z = 1) follows from o" = a. If o" = 0, 
we have Z = Z. Symmetry II in connection with the trivial renormalization 
of c shows that / is renormalized with the same Z factor as r: Zy — Z"^. 

The RGE (43) now contains four Gell-Mann-Low functions j3i corresponding 
to the four dimensionless parameters v, w, and p. A one- loop calculation 
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gives 



3u 2(5 + 5p + 2p2)t; (7 + 8p + 3p2)w 



+ 



+ 



2w 



.1 + p p{l + p) p{l + p) 
I3y = (-£ + 2«;)v 



+ 0(2-loop) 



S'U 4?; 

^^T" (l + p)3 



/9„ 



— £ + — 



2(l + p)3 
2(3 + 2p + p2)v (5 + 5p + 2p2)w 



;i + p)3 + (i + p)3 

2v (7 + Ip + i)^)w 



+ 0(2-loop) 



+ 



8 (l + p)3 4(l + p)3 



0(2-loop) 



P- 



(173a) 
(173b) 

(173c) 
(173d) 



The RG flow functions leading to the anomalous dimensions of the flelds s 
and s are found to be 



7= S 



7 + 7 



7 if (T = 

if (7 7^0 

4t; 



(3 + p)w 



u 

"2 + (l + p)2 (l + p)2 



0(2-loop) 



(174a) 

(174b) 



At a non-trivial IR-stable fixed point with all /?j = and and p* both 
different from and oo, one finds from Eqs. (173b) and (173d) the exact 
statements = e/2 and C* = 0, and Eq. (47) yields the asymptotic scaling 
properties of the Green's functions, 



N,N;M,M\ 



(175) 



Consequently, the DP-C universality is characterized by the exact critical ex- 
ponents 

In Eq. (175) we have 



5g 



(M + M) + {Np + 7V/5')J d/2 , 
/3' if (7 = 
1 if (7^0 ' 
/3'=(d + 77')/rf, V = 7*- 



= < 



(176) 

(177a) 
(177b) 
(177c) 



Thus, there is merely one unknown independent critical exponent tj' in the 
DP-C universality class. 
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To order e, one finds the following IR fixed points as zeros of the Gell-Mann- 
Low functions 

1 27 

cr = : P* = - , u^ = 2e , = — e , = ; (178a) 

2 o4 

(7>0: p, ^ (2 + Vsf^^ + (2 - Vsf^^ - 2 , (178b) 

4 1 + P* 1 - 5p* 

= -r e, ^ — - — e, = £; (178c) 

2 + 4 

cr < : run-away flow. (178d) 

The run-away solution in the case a < eventually violates the stabihty 

bound u > 2w in the course of the flow, which indicates the emergence of 
a first-order transition is this situation. In the other cases, we find for the 
remaining exponents to one-loop order: 

(7 = 0: 77' = -e/8 , (179a) 
a>0: 77' = — . (179b) 

All scaling laws now follow from the scaling of the Green's functions (175). 
As important examples we note the initial slip exponent 6i — —v'f^ 
9i — —r]'/4: for cr = and cr > 0, respectively. 



5.5 Directed Percolation and Quenched Disorder 



According to the equation of motion (162a), the DP-C processes can be un- 
derstood as models of systems where the critical control parameter r itself 
becomes a dynamical variable r + /c(r, t) subject to a diffusive process. For 
the case that a DP process evolves in a disordered medium, the field c repre- 
sents a static but spatially inhomogeneous random quantity with short-range 
Gaussian correlations c(r) c(r') oc 6{r — r'). It is easily seen that the influence 
of the disorder on the other parameters of the model and beyond-Gaussian 
correlations are IR-irrelevant. In contrast to equilibrium systems, where the 
normalization factor of the probability distribution, namely the partition func- 
tion, is itself a functional of the disorder, the description of dynamics by means 
of an exponential weight exp(— ^c) with a response functional J^c that depends 
on the disorder fleld c(r) does not require a normalization factor. Hence, aver- 
aging over c can be easily performed directly on the statistical weight directly, 
exp(— jTc) =: exp(— JT"), without invoking any additional procedures such as, 
e.g., replication. One thereby obtains the following effective action for the 
evaluation of path integrals averaged over the randomness [73]: 



J = ld'^r[ldt~s[d, + A(r - V^) + ^{s -s)]s-^f[J^ 



dt ss 



(180) 
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Prom here, the calculation of the renormalizations and the RG flow functions 
proceeds in the usual manner. To one-loop order, one finds the Gell-Mann-Low 
functions 

(-£ + Y-10^;)ii, p,^ (^-e + ^-8vy, (181) 

where v is the renormalized dimensionless coupling corresponding to /. It is 
now easily seen that the flow equations (181) allow only for runaway solutions 
as / — >^ in the physical region u > 0, v > 0. The fixed point of the pure system 

— 2s/3, — is unstable, as is the non-physical fixed point u^, — —4s/9, 

— —1/6. There is a stable fixed point, namely — 0, — —1/8, but it 
is in the non-physical region as well. The runaway RG trajectories render the 
perturbation expansion useless and, perhaps, indicate a more complicated crit- 
ical behavior than just simple power laws. Indeed, the simulations by Moreira 
and Dickman [74] show logarithmic critical spreading instead of power laws, 
and seem to display violation of simple scaling [75] . Apparently the continuum 
field-theoretic treatment is not capable of capturing the strong disorder limit 
of DP, which may be dominated by rare localized excitations. A promising 
alternative approach has recently been developed that employs a real-space 
RG framework specifically tailored to strongly disordered systems [76] . 
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